APPLICATION OF CURRENT ALGEBRA OR CHIRAL SYMMETRY 

TO TAU HADRONIC DECAYS* 



1. Introduction 



L. BELDJOUDI 

Centre de Physique Theorique, Ecole Polytechnique 
91128 Palaiseau, France 
E-mail: beljoudi@orphee.polytechmque.fr 

ly-s ' and 

Q^, T. N. TRUONG 

' Centre de Physique Theorique, Ecole Polytechnique 

O ' E-mail: truong@orphee.polytechnique.fr 

CD ■ 

Q: 

■ Abstract 

I r — > vrvrz^, r — > t:Ku, t Krjv, r Snu and r — > ttttKi^ have been 

' investigated using chiral symmetry with dispersion relation in agreement with 

Tlj" ■ the unitarity condition. 
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-p^ ' Current Algebra was invented a long time ago to study phenomena 

p • involving emission of soft pions (and kaons) and also the chiral symmetry 

p ■ breaking effects. It was soon discovered that the current algebra low en- 

D . ergy theorems (LET) can be obtained by the effective lagrangian method 

', at the tree graph order. Using it, we can calculate matrix elements of 

.!-H i physical processes by the Feynman perturbation method in a straight- 

^ I forward way. Because of the chiral properties of the Nambu-Goldstone 

' bosons, the effective lagrangian involves the derivatives of the pion fields 

which make the theory non-renormalisable. Chiral Perturbation Theory 
(xPT) is a well defined perturbation procedure which can take into ac- 
count in a systematic way of higher orders (higher loops) and of the chiral 
symmetry breaking effects. At any order, these relations form low energy 
theorems of QCD. Despite the beauty of the method, it has limitations in 
phenomenological applications. We distinguish three main problems with 
this theory: 

1) Because of the non-renormalisability of the theory, the number of 
parameters increases as we calculate higher loops. 

2) xPT is effectively a power series expansion in terms of pion mo- 
menta, hence it cannot take into account of the resonance effects (Breit 
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Wigner form cannot be expanded in a convergent power series of mo- 
menta). It can only explain the low energy tail of the resonance. 

3) Because xPT is a perturbation theory, it only satisfies the unitarity 
relation order by order. This approximation is not good enough for strong 
interaction or for resonance physics. 

Because of these problems, at first sight, xPT cannot be used to study 
Tau physics where most of the hadronic modes are dominated by reso- 
nances and where the pions emitted are not soft. The remedy for these 
problems were given a long time ago even before xPT became popular. 
One method is to use the Current Algebra soft pion theorems and use dis- 
persion relation to take into account of the unitarity. The second method 
is to extend the validity of the xPT by resumming the perturbation series 
in order to take into account of the unitarity relation. The third method 
is to assume the vector meson dominance in the elcctroweak processes 
and requiring at low energy or at zero momentum transfer the low enegy 
theorem is recovered. In fact, the three methods are equivalent as they 
all satisfy the LET and the unitarity relation to a good approximation. 
Using the basic Equal Time Commutator Relations (ETCR) as postulated 
by Gell-Mann, 



= ^r"'A^,{x) (3) 

and are the vector and axial currents generated by the approximate 
symmetry SU{S)l x 5'C/(3)r of QCD, Qy and are the corresponding 
charges. 

Together with the basic Current Algebra formula for a soft pion emis- 
sion : 

]im {7r'^{k)B\0mA) = -f [g^,O(0)] \A) (4) 

Eq(4) sets the scale for the relevant matrix clement. In the disper- 
sion relation language it can be used as a substraction constant for a 
substracted dispersion relations. 

More explicitly, let us denote the Current Algebra result as Alet{0), 
thus the substracted dispersion relation can be written as follows: 

A{s) = A,et{0) + - Jdz /"^ ^^^l . (5) 
TT J ziz — s — le) 

so ^ ' 

Because of the substracted dispersion relation, the high energy contri- 
bution in the dispersion integral, which is difficult to calculate reliably, is 



(1) 



(2) 



suppressed. This result is in contrast with the physics done in the 60's 
(e.g. Bootstrap) which involves unsubstracted dispersion relation with un- 
controllable approximation (because the high energy contribution is not 
suppressed) . 

In the physical world, where the pion are not soft, corrections must 
be made. First the chiral symmetry breaking effect can be taken into 
account approximately by hand. The second correction to be made is to 
impose unitarity relation which is crucial for hadronic processes. This can 
be done either by strong interaction dynamics models e.g. Vector Meson 
Dominance or some unitarisation schemes. 

Some of these relations were calculated in 19780. It is desirable to 
update these calculations with a more sophisticated technique in order 
to compare with more recent experimental data and also to the xPT 
technique. As an illustration we present applications to r ^ tttiu, t —>■ 
ttKu, t —>■ Krju T —>■ npu, r —>■ ttK*v and r Kpv decays. 

2. r — s> TTTTZ/ decay 

Using the CVC hypothesis and Lorentz invariance, one can write straight- 
forwardly the following matrix elements: 



At zero momentum transfer we have the following normalization: /(O) = 1 
and r(0) = m^. Assuming the elastic unitarity condition, we deduce the 
following relations: 



Hence f{s) must have the phase of the P wave tttt phase shift, and 
r(s) the S wave phase tttt shift. 

The general solutions to this equation are well known, they are of the 
Muskhelishvili Omnes typeH: 
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where 
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^^o(^) = eM- [^^^) (12) 
TT J z[z — s — le 

Pf and Pr are polynomials normalized to unity at s = which de- 
termine the high energy behavior of the form factors. They could also 
represent the low energy contribution of the higher mass intermediate 
states to the form factors. In the following we assume the dominance of 
the elastic unitarity relation and hence we set Pf{s) = Py-{s) = 1. 

At one loop chiral perturbation theory we have: 
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f{s) =1 + f s - {^^Tif'r [s - Ami){h{s) - h{0)) + 4mih'{0)s (13) 
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If we make a bubble summation consistently with the elastic unitarity 
condition, we obtain: 



where 



l-'fs + (967r/2)-i [(s - Aml){h{s) - h{0)) + Amlh'{0)s] 

m 



l-!fs + (16vr/2)-i (s - ^^m-s) - hiO)) + ^^sh'iO) 



(15) 



his) ^ lpis)H^±^^) (16) 

p{s) is the phase space factor. 

Eq(13) can also be obtained by the diagonal Fade approximant method 
by using Eq(ll) and (12). This approximation is identical to the Vector 
Meson Dominance approach of Gell-Mann-Sharp and Wagner where the 
self energy correction is taken into account to make the vector meson 
instable and with pmi coupling given by KSRF relationS. 

An alternative method to calculate the scalar and vector form factors 
is to use the S and P wave tttt phase shifts and the Omnes representation 
(see Eq.(9) and Eq.(lO)). For more details see Ref.i. The P wave rni 
phase shift is given in Fig.(l). 

Because there are no experimental information on the scalar form fac- 
tor, we only compare the vector form factor with the experimental data. It 
is seen that the agreement between theory and experimental data is satis- 
factory although the peak values of the experimental form factor squared 
at the p mass is about 40 as compared with the theoretical calculation 
value 32 or an error of the order of 20% as it can be seen in Fig. (2). 
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Figure 1: The solid line I — 1, I — 1 pion scattering phase shift calculated from 
unitarized ChPT Dashed line corresponds to the similar phase shift where the left 
hand cut is neglected. 
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Figure 2: The solid line is the vector pion form factor squared calculated from the 
Omnes representation. The dashed line represents the same quantity calculated by the 
Fade Method 



This discrepancy is probably due to the inelastic effect of the chan- 
nel as was previously pointed outi. 

We can also calculate the vector and scalar pion rms radii using the 
following formula: 

6 5ldz 
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Numerical integration gives (r^) = 0.40/m^, and (r|) = ^Aljvn? com- 
pared to the experimental value (r^) = 0.439 ± 0.03/m^, and the (r|) = 
O.bfm^ which is obtained from the experimental irK scalar radius and 
from SU(3) symmetry. The agreement between experimental data and 
theoretical calculations is satisfactory. 



3. T ^ ttKu decay 

The most general r ttKu decay amplitudes are given in terms of 
two form factors: 

{Tr^K-\V;-^'{0)\0) = Ms){p2-pi), + f2{s){pi+P2), (17) 

where pi, and p2 are, respectively, the pion and kaon momenta, and s = 
[Pi + ^2)^ is the time-like momentum transfer and is the vector 

current operator with the superscript indices referring to the SU(3) octet 
currents. fi{s) is the P wave irK form factor, /2(s) is a linear combination 
of S and P states as can be seen by taking the divergence of Eq (15): 

gis) = -^{n'K-\^X-''im) = « " OMs) + sU-s) (18) 

g(s) is therefore a pure scalar which describes the S wave nK form factor. 
g(s) measures the SU (3) violating effect because, in the exact SU(3) limit, 
the vector current is conserved. We expect therefore in the r -kKv 
decay, the P wave form factor fi{s) dominates. 

Because of the octet current hypothesis the two channels Tr^K~ and 
n'K^ matrix elements are related by the Clebsh Gordon coefficient 

(7r-ir°|V;^"^5(0)|0) = V2(7r°ir-|V;^-*^(0)|0) (19) 

Using the Ademollo-Gatto theorem /i(0) = and hence 

g{0) = {mj^ — m^)/-\/2 for n^K~ system. 



Using the standard current algebra technique and the SU{2)lX SU{2)r 
commutation relation by taking the pion momentum pi soft we have the 
well known Callan-Treiman relation^: 

h{m\) + h{ml) = (20) 

where fx and /tt are, respectively, the K and tt decay constants 
fKlh = 1-22. 

By evaluating Eq(16) at t = m\ and noting that f2{jnK) is propor- 
tional to ml/mj^ we have: g{m\) ^ g{0)fK/ fn 

The elastic unitarity condition, which should be valid in the physical 
region of the r — >• ttKu decay gives: 

/m/i(s) = his) exp -zSy^ sin (21) 
Img{s) = g{s) exp —16^^ sin (22) 

where and are respectively the phase of S and P wave 1=1/2 nK 
scattering amplitude. The general solutions are given by the Muskhel- 
ishvili Omnes representationcl: 

+00 j;1/2j 

/(.) = /(0)exp(^ / ' ' . ) (23) 
TT J z[z — s — le 

g{s)=g{^)^M- \ , ' . ) (24) 
TT J z[z — s — %e 



At one loop chiral perturbation theory we have: 

/(s) = /(O) + t'^^s) (25) 
g{s) = g{^)^g'°°^{s) (26) 

/'°°P(s) and g^°°^{^s) are given elsewhereB. If we make the usual bubble 
summation, consistent with the elastic unitarity condition, we obtain: 



/'«>P(3)//(0) 



From the expression for fi{s), the phase of the form factor which is 
identical to the P wave phase shifts of irK scattering amplitude, can be 
calculated using the experimental value of (r^) = 0.34 ± O.OSfm^. Using 
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Figure 3: The calculated P wave ttK form factors (solid/dashed/dot-dashed curves) 
corresponding respectively to, the Omnes representation using the ttK phase shift 
calculated from the unitarized one loop CPTh, the [1,1] Fade approximant, CPTh 
prediction 



this value we have mx* = 810 ± 30 MeV which agrees with the exper- 
imental data niK* = 892MeV. Its width satisfies the following modified 
KSRF relation! 



.2 



Using the experimental value ttik* = 892MeV the numerical result of the 
right hand side of Eq(26) is 55 MeV, compared to the experimental value 
of 49.8 ± 0.8 MeV. In Fig. (3), we give a graphic representation of 
squared. 

The branching ratio B.R = jg 1.0% and is in agreement with 

the experimental result of B.Rexp. = (1-4 ± 0.2)% 

Because the S wave nK scattering length does not vanish, g(s) has a 
square root threshold singularity at the threshold (the derivative of g(s) 
is discontinuous at this point) as it can be seen in Fig. (4). 

The scalar form factor contributes very little to the ttK spectrum owing 
to the fact that it appears as a square of the amplitude. The forward 
backward asymmetry, being proportional to the amplitude, is reasonably 
large. It is about 10% in the K* resonance region where the number of 
events is maximum as it can be seen in Fig. (5). 




S [GeV'^2] 

Figure 4: The calculated S wave ttK form factors (solid/dashed/dot-dashed curves) 
corresponding respectively to, the Omnes representation using the TtK phase shift 
calculated from unitarized one loop CPTh, the [1,1] Fade approximant , CPTh pre- 
diction. 

The forward-backward asymmetry could be a useful quantity for study- 
ing the relative phases of the S and P waves. 

The Omnes representation with the S and P wave J = 1/2 ttK scatter- 
ing amplitude studied in Ref.0 (and given in Fig. (6)) yields a branching 
ratio of 1.15% for r ttKv decay, in a better agreement with the exper- 
imental value 1.4 ± 0.2%. 

4. r ^ Krjv 

In this process elastic unitarity condition does not apply. One should 
consider the two channels TtK and Krj in the final state interaction. Chiral 
Perturbation Theory cannot describe the Krj scattering even at thresh- 
old since the K* resonance occurs below Krj threshold. According to 
the equivalence between the Pade Method and the VMD approach (as 
it was pointed out in the pion form factor calculation), we will use in 
the following the K* dominance for studying the r — > Krjv decay. VMD 
method introduces the inelastic effects through the K* self energy. Using 
the exprimental value of r — >■ K*!^ branching ratio and SU(3) symmetry, 
we predict a branching ratio of 1.6 x 10^^. The experimental value of 



S [GeV'^2] 



Figure 5: Prediction for the forward-backward asymmetry ApB defined in eq(5) as 
a function of the ttK invariant mass squared. The dashed/solid curves correspond 
respectively to the calculations with/without the left hand cut ofirK scattering ampli- 
tude. 
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Figure 6: The solid line represents the 1=1/2, 1=1 nK scattering phase shift calculated 
from the unitarized CPTh. The dashed line corresponds to the similar phase shift 
when the left hand cut is neglected. The dot-dashed line is the CPTh prediction phase 
(which is not the same as the phase shift due to the violation of the full elastic unitarity 
relation in this method). 



r —>■ Krjv branching ratio is (2.5 ± 0.5) x 10 ^. 



5. r ^ upv, T — s> ttK*^ and r Kpu 

The most general matrix element can be written as: 

{7,-{k)p\p)\Al-^\m = h{Q')e,+e.k [{k + p)MQ') + {k-p),h{Q') 

(29) 

where = {k+pY and e is the polarisation vector of p. fi, /2, and /s are 
complex form factors and are only functions of Q^. Current algebra soft 
pion theorem, which is obtained by taking the limit A;^ — 0, gives only 
information on /i but not on the other two form factors. In an explicit 
model, it was shown that they contribute little to the r —>■ irpu. Interested 
readers are refered to the original article 0. (We assume here that the decay 
constant of vr' is sufficiently small and hence can be neglected). Using the 
standard low energy current algebra theorem and taking the limit A;^ ^ 
we have: 

\im{n-{k)p'{p)\Al-'^\0m = -V2^e,{p) (30) 

where = 93MeV, and fp is defined by the rate of p ^ e~^e~. Using 
the experimental data we obtain, fp = O.llSGeV^^. This value of fp is 
equivalent to writing approximately the pion form factor as 
Ft,{s) = m'^p{l + 6s/m?p)/ {ni^p — s — i'mpTp{s)^. A good fit to the experi- 
mental data is obtained with 5 = 0.2. In fact, the more general form of 
Eq(3) reads 



hmJvr-(A;)vr+(gi)^-(g2)l4-^'(0)|0) = -^F.(.)(gi - q,)^ (31) 

In 



For convenience we shall first use Eq(3). The 3tt matrix element below 
the pir threshold can be straightforwardly obtained from Eq(28). Using 
Eq(26) in (27) we have: 

f,{ml) = -V2^ (32) 

In 

Let us start with the narrow width approximation for the Ai propagator. 
Using Ai dominance for the form factor we have: 



/i(g^) = -V2if ^^^_2 (33) 



-Jp im\-ml 
' fn m\ - g2 

The generalisation of Eq[30] to take into account of the unstable nature 
of Ax can be straightforwardly made. Using the A^ dominance hypothesis 



for the axial current, the general expression for /i(Q^) is: 

f mi — mi — IT (mi) 
= -v^^^^^f ^ 34 

where we use the standard prescription for describing an unstable particle, 
with 77 (Q^) being the self energy operator of the Ai resonance and is 
obtained by the bubble summation of the intermediate states, similar to 
the treatment of the W and Z propagators in the standard model. In order 
to have the usual Breit Wigner description of a resonance, we must make 
a twice substracted dispersion relation with Re[7i{m\)] = Re[7i'{m\)] = 
where rriA is the Ai mass: 



(Q^ - m\)^ 7 Im['K{z)]-Im[K{m\)]-{z - m\) Im[7r' {m\)] 



Re[rr{Q^)] J^ J"^' pjdz 



where we define the 7r°p+Af vertex as gAp-K^{A).e{p), A [Q^., m^, rn^ = 
(Q^ — (mp + m^)^)((5^ — {mp — m^)^), and P stands for the principal part 
integration. 

If we assume that the acceptance correction to the ALEPH data was 
negligible, our best fit to the Svr spectrum gives the following ranges of 
rriA and F^: = 1.24 ± 0.02GeV, Ta = 0.43 ± 0.02 GeV. 

Using these results, the calculated branching ratio for r — > Svri/ is 
19 ± 3%. The central value corresponds to our best fit which is shown in 
Fig. (7). This value agrees with the ALEPH branching ratio 19.14±0.48± 
0.44%. 

A similar treatment for r 'kK*v and r — > Kpv decays can be done 
using the axial vector dominance of Qi(1270) and (52(1400) (for more 
details see Ref.i). Numerical calculation gives: 

BR{t ^ p^K-v) = 0.1% 
BR{t n-K°*u) = 0.4% (36) 

These results are in good agreement with the values of TPC/Two-Gamma 
collaboration: B {r ^ K*^ti~v^ neutrals) = 0.51 ± 0.2 ± 0.13 and 
B {t ^ K^Ti^Ti^u) = 0.7 ± 0.2. The K^p^u mode is therefore consistent 
with zero. 



6. Conclusions 



To conclude we would like to outline some main points: 
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Figure 7: Our best fit for the r Siriy spectrum corresponding to rriA = 1-24 GeV, 
= 0.43 GeV 

(i) Perturbation theory satisfies unitarity order by order and that is 
known as perturbative unitarity. 

This method is sufficient to describe physical phenomena if the cou- 
pling constant is small like in QED. But it fails if there is a bound state 
or a resonance state: in QED the positronium bound state is solved by 
Bethe-Salpeter equations (Ladder type summation). 

(a)- Notice that in the Standard Model, the unitarity is preserved by 
the bubble summation of the propagator. 

-In the Nambu-Jona-Lasinio Model the bubble summation is done to 
satisfy the unitarity. The Higgs and Nambu Goldstones are generated as 
bound states due to four fermion interaction. 

(in) Low energy Nucleon-Nucleon interaction has resonance in the sin- 
glet S state and a bound state (Deuteron) in the triplet state. This prob- 
lem was solved by H. Bethe in the 50's using effective range expansion 
which preserves the unitarity of the S matrix. More explicitly, instead of 
expanding the amplitude like in Chiral Perturbation Theory, the Effective 
Range Theory makes an expansion of kcoi{5) 

A;cot(5) = - + -roA;2 + ... (37) 
(X 2 

with a and ro being respectively the scattering length and the effective 
range. 



This expansion, unlike xPT, preserves unitarity. 

(iiii) In this paper wc have shown that Chiral Perturbation Theory 
should be modified to describe the main problems of r physics. We have 
extended the range of validity of Chiral Perturbation Theory using the 
Pade Approximant Method or VMD hypothesis. Both methods satisfy 
unitarity. A good agreement with the experimental data is obtained for 
the exclusive hadronic modes of r decay. 
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